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ABSTRACT 

Two different definitions of rotational relaxation times were used to calculate rota- 
tional collision numbers for pure polar gases. One definition was based on the relaxation 
equation for rotational energy while the other was based on the connection between relax- 
ation time and volume viscosity. For the purposes of the calculation, a polar molecule 
was taken to be a point dipole imbedded in a hard core and the interaction between mole- 
cules was confined to a plane. A classical calculation showed that the two definitions are 
equivalent, at least through the third order of a perturbation calculation. Both defini- 
tions predicted that increasing the temperature and decreasing the dipole moment and 
moment of inertia should increase the rotational collision number. The results are in 
qualitative agreement with experiment insofar as the dipole moment and moment of inertia 
dependence are concerned. No conclusion can be made concerning the correctness of 
the theoretical temperature dependence since the experimental temperature dependence 
is somewhat in doubt. 
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A COMPARISON OF ROTATIONAL COLLISION NUMBERS 
OBTAINED FROM DIFFERENT DEFINITIONS 
by Frank J. Zeleznik 
Lewis Research Center 

SUMMARY 

Two different definitions of rotational relaxation times were used to calculate rota- 
tional collision numbers for pure polar gases. One definition was based on the relaxa- 
tion equation for rotational energy while the other was based on the connection between 
relaxation time and volume viscosity. For the purposes of the calculation, a polar 
molecule was taken to be a point dipole imbedded in a hard core and the interaction be- 
tween molecules was confined to a plane. A classical calculation showed that the two 
definitions are equivalent, at least through the third order of a perturbation calculation. 
Both definitions predicted that increasing the temperature and decreasing the dipole 
moment and moment of inertia should increase the rotational collision number. The 
results are in qualitative agreement with experiment insofar as the dipole moment and 
moment of inertia dependence are concerned. No conclusion can be made concerning the 
correctness of the theoretical temperature dependence since the experimental tempera- 
ture dependence is somewhat in doubt. 


INTRODUCTION 

In comparison with vibrational relaxation, relatively little experimental data on ro- 
tational relaxation exist in the literature. This is particularly true of rotational relaxa- 
tion in polar gases, and, until quite recently, experimental values have been almost 
nonexistent. However, recent experimental studies have yielded rotational relaxation 
times, or equivalently, rotational collision numbers for a number of the more common 
polar gases (refs. 1 to 10). Some of the available data is summarized in table I. The 
data (ref. 9) for the molecule NO have not been included because it is considerably less 
polar than the molecules of table I, and hence the dipole -dipole potential is probably not 
the dominant interaction. In addition, a z jjqt for ammonia, determined by Petralia 
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(ref. 11), has not been included since it is approximately four times larger than the cor- 
responding value in table I. In a private communication, Dr. R. S. Brokaw of this 
laboratory expressed the opinion that the Z r q t for HgS that is given in table I appears 
unrealistically large when compared with similar molecules. However, I have included 
it in the tabulation because it is the only data available for this molecule. Except for 
the fluoromethanes, the experimental collision numbers have been obtained from an 
analysis of either acoustic absorption or thermal conductivity measurements. The 
fluoromethane data were derived from nuclear magnetic resonance spectra (ref. 10). 

It is by no means certain that the three experimental techniques give comparable values 
for the rotational collision number The two molecules HC1 and HgO can be used 

to compare the results from thermal conductivity and acoustic measurements. These 
two methods appear to give consistent results as is shown in figure 1. This agreement 
is quite important since thermal conductivity measurements coupled with the theory of 
Mason and Monchick (ref. 12) promises to be the most convenient way of obtaining Z r q T> 
The NMR experiments also appear to give results that are comparable to those from 
thermal conductivity measurements, at least for spherically symmetric molecules 
(ref. 10). 

There are three obvious parameters that could affect the rotational collision number 
of a polar gas: (1) the magnitude of the permanent electric dipole moment, (2) the mass 

TABLE I. - PROPERTIES OF POLAR GASES 


Molecule Electric Molecular 


ch 2 f 2 


dipole 

moment, 

2 

(esu)(cm ) 

1.081X10" 1£ 

1.085 

1.83 

1.837 

1.844 

1.861 

1.477 

1.509 

1.477 


weight, 


36.469 
37. 476 
20. 008 
21. 014 
18. 016 
20.029 

17.032 
20.050 

18. 033 


(a) Molecular properties 

o 

Moments of inertia, (g-cm ) 


1. 0224x10" 


1.833 
4.4140 
8. 7960 
4.4140 


64.063 95.351 

70.019 81.08 

52. 026 17. 1 

34. 034 5. 52 

34.080 5.8008 


2. 6431x10" 

5. 1404 

1.33563 

2. 54140 

1.9180 

3.841 

2.8087 

5.4129 

2.8159 

81. 322 
81. 08 
79. 1 
32.854 


2.9404x10" 

5.674 

2.8087 

5.4129 

2.8159 


13.806 
148.67 
90. 7 
32.854 
2 . 7040 


Average 
moment of 


2. 6431X10" 
5. 1404 

1. 33563 

2. 54140 
1.9603 
3. 7827 
3.3438 
6. 5406 
3.3486 

63.493 
103.61 
62. 3 
23.74 
3.8693 


Moment of 
inertia to 
molecular 
weight ratio, 
I/M 


0. 725x10" 

1. 372 


.991 
14. 797 
11.975 
6.975 
1. 135 
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TABLE I. - Concluded. PROPERTIES OF POLAR GASES 
(b) Experimental rotational collision numbers 


Molecule Tempera- 
ture, 


Rotational Data 
collision from 
number, refer - 


273. 

2 

7. 

0 

300, 

1 

6. 

2 

328. 

5 

4. 

6 

374. 

8 

3. 

6 

423. 

1 

3. 

2 

471. 

4 

3. 

0 

300. 

1 

2. 

6 

328. 

5 

2. 

3 

374. 

8 

2. 

0 

423. 

1 

1 . 

9 

471. 

4 

1 . 

9 

373. 

, 75 

9. 

5 

422. 

30 

9. 

1 

373. 

75 

4. 

1 

422. 

30 

3. 

8 

323. 

2 

4. 

0 

381. 

2 

2. 

7 

426. 

1 

2. 

5 

478. 

0 

2. 

3 

525. 

6 

2. 

0 

381. 

2 

1 . 

6 

426. 

1 

1 . 

5 

478. 

0 

1 . 

3 

525. 

6 

1 . 

2 

300. 

0 

2. 

3 

329. 

1 

2. 

1 

374. 

6 

1 . 

9 

424. 

2 

1 . 

7 

474. 

5 

1 . 

6 

311. 

2 

1 . 

665 

333. 

2 

1 . 

700 

353. 

2 

1 . 

709 

373. 

2 

1 . 

866 

393. 

2 

2. 

000 

413. 

2 

2. 

073 

433. 

2 

2. 

144 

453. 

2 

2. 

299 

473. 

2 

2. 

628 


Method Molecule Tempera- 
ture, 


Rotational Data 
collision, from 
number, refer - 
Zrot ence 


CHF 3 

300 

2.01 

10 

N 

CH 2 F 2 

300 

1. 54 

10 

N 

CH 3 F 

300 

3. 61 

10 

N 

h 2 s 

298. 2 

31. 2 

6 

A 


T, thermal conductivity; A, acoustic absorption; N, nuclear magnetic resonance. 
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Figure 1. - Comparison of rotational collision number as obtained 
from acoustic absorption and thermal conductivity measurements. 


distribution as characterized by the ratio of an arithmetic mean moment of inertia I 

to the mass M, and (3) the temperature T. The qualitative effect of these parameters 

can be deduced from the values of and the molecular properties that are given in 

table I. Thus, the pair of molecules HgO and HgS have approximately the same value of 

I/M but appreciably different dipole moments; a similar relation exists between the 
15 

pair DgO and NHg. The data for these two pairs of molecules indicate that an increase 
in the dipole moment produces a decrease in Z^Q,p. This effect was first noted by 
Brezeale and Kneser (ref. 8) when they compared the Z^Q,p for HC1 with that of sim- 
ilar nonpolar molecules. To determine the role of mass distribution in the rotational 
relaxation process, we can compare the pair of molecules HC1 and DC1 since they have 
very nearly the same dipole moments. This is also true of the pair H^O and DgO and 
the set of molecules NHg, ^NHg, and NDg. These comparisons indicate that an increase 
in I/M produces a decrease in Z^q^* Finally, the data of table I can be examined to 
determine the temperature dependence of Z^q^. All the molecules, with the exception 
of SOg and NHg, have a that decreases with temperature. The Z r q t for SOg 

increases with temperature. For NHg, the experiments of Baker and Brokaw (ref. 2) 
gave a that decreases with temperature, but the similar experiments of 

Srivastava and Das Gupta (ref. 5) showed the opposite results. It is interesting to note 
that only a small change (2 to 3 percent) in the thermal conductivity values of Srivastava 
and Das Gupta is sufficient to make Z R q t a decreasing function of temperature (ref. 5). 
This sensitivity makes the temperature dependence of Z r q t uncertain. 

The experimentally deduced effects of the parameters I/M and T cannot be com- 
pared with theoretical predictions since no theoretical calculations of Z r q t for polar 
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gases exist. However, some classical calculations have been made for a number of 
nonpolar gas models. Sather and Dahler (refs. 13 and 14) have examined the rough 
sphere, the spherocylinder, and the loaded sphere gas molecules. For the rough sphere 
the find 


2 

Z ROT “ ( 4i /Mct 2 ) /(l + 41 /Mct 2 ) 

where I is the moment of the inertia, M the molecular mass, and a the molecular 
diameter. Because 4I/Ma « 1, the decreases with increasing I/M. For 

the spherocylinder and the loaded sphere, their calculations show the opposite depend- 
ence, that is, Zjy-j.p increases with increase in I/M. Widom’s calculations (ref. 15) 
for rough sphere molecules being relaxed by an inert gas give a dependence on I/M 
similar to Sather and Dahler’ s rough sphere results. Using a somewhat more realistic 
potential but restricting collisions to a plane, Parker (ref. 16) showed that the rotational 
collision number increased linearly with I/M. Finally, Brout’s (refs. 17 and 18) quan- 
tum mechanical calculations imply a decreasing Z^q^ with increasing I/M. The 
theoretical temperature dependence of Z^Qrp, like its mass distribution dependence, is 
a strong function of the model used in the calculation. For rough spheres, spherocylin- 
ders and loaded spheres, Sather and Dahler (refs. 13 and 14) find Z^q^ to be inde- 
pendent of temperature. However, when the rough spheres and the spherocylinders are 
surrounded by a square well potential, the rotational collision number becomes a rapidly 
increasing function of temperature. Parker (ref. 16) also obtained a Z^q^ that in- 
creases with temperature although at a somewhat greater rate than the square -well 
models. 

Recently Nyeland (ref. 19) reapplied Parker’s model for the calculation of Zj^q^, 
but used a different definition for Zj^q^; the two possible definitions are discussed in 
the next section. In contrast to Parker’s results, Nyeland's Zj^q^, decreased with an 
increase in the moment of inertia and had a smaller rate of increase with temperature. 
This difference might be due either to the use of different definitions or to different 
approximations in the calculation. 

In view of the present uncertainties concerning the rotational relaxation of polar 
gases, it seemed desirable to attempt a theoretical treatment of the problem. It would 
be preferable to carry out a quantum mechanical treatment, and, indeed, Cross and 
Gordon (ref. 20) have recently evaluated the first Born approximation to the scattering 
matrix for scattering from point dipoles. However, they make no attempt to account 
for short range interactions since their primary concern is with the effect of resonant 
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scattering on the total scattering cross section. It has been conjectured by Mason and 
Monchick (ref. 12) that resonant scattering is also a significant effect in rotational re- 
laxation. However, based on experimental results, Baker and Brokaw (ref. 1) suggest 
that the effect may be classical rather than quantum mechanical. For this reason, it 
would be desirable to see what extent classical mechanics can account for the experi- 
mental observations. An additional objective is to ascertain the effect of using different 
definitions of relaxation times. 

The calculation I shall describe is essentially a classical perturbation calculation of 
the change in rotational energy produced by the scattering of permanent electric dipoles. 
The molecules are constrained to a plane and no out-of -plane forces are permitted. 

This sacrifice of some of the physical reality is necessary in order to make the calcu- 
lation mathematically tractable. The calculation is by no means trivial even with this 
assumption. The interaction potential and the details of the computation differ consider- 
ably from Parker's, although both are two dimensional. I shall employ an interaction 
picture of classical mechanics as described by Garrido (refs. 21 and 22). For con- 
venience, pertinent features of Garrido' s method are summarized in appendix A. 


ROTATIONAL RELAXATION TIMES 

Two definitions of relaxation times have appeared in the literature. These two 
definitions have not been shown to be equal; however, they are often used interchange- 
ably. The first of these definitions is a phenomenological one and corresponds to the ad 
hoc introduction of an additional equation (and an adjustable parameter) to supplement 
the usual hydrodynamic -thermodynamic equations of continuum mechanics. This sup- 
plementary equation was initially introduced to explain acoustic absorption experiments. 
The new parameter, known as a relaxation time, is defined by the supplementary 
equation. If we restrict our attention to rotational energy, then the rotational relaxa- 
tion time r R is defined by 


dE R _ E R< T > - E R< t > (1) 

dt t R 

In this equation E-^ is the contribution of rotational energy to the thermodynamic inter- 
nal energy. The equilibrium value of the rotational internal energy is Ej^(T) while the 
instantaneous, nonequilibrium value is E^(t). 

The alternate definition of a relaxation time is based on the fact that at low acoustic 
frequencies it is impossible to experimentally separate the absorption of acoustic energy 
due to viscous -thermal effects from that due to a relaxation phenomena (ref. 23). The 
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absorption can be corrected for the effects of shear viscosity and thermal conductivity 
using experimentally determined values of shear viscosity and thermal conductivity. 

It cannot, however, be corrected for the absorption due to volume (or bulk) viscosity 
since there is no known experiment that independently measures this transport coeffi- 
cient. Consequently, little is known about the volume viscosity and it has become com- 
mon practice to equate the volume viscosity absorption with the relaxation absorptions. 
This leads to a relaxation time, r R , defined in terms of the volume viscosity k by the 
equation 



where v is the number density of particles, C y is the heat capacity at constant volume, 

and C.. is the contribution of internal motions to C y . The contribution to is 

only from rotational motion in our case. Equation (2) can be obtained by combining 

equations (38) and (74) of the paper by Wang Chang and Uhlenbeck (ref. 24). Wang Chang 

and Uhlenbeck have deduced an expression for k by generating the normal solution of 

the Boltzmann equation and I shall use the classical equivalent of their semiclassical 

expression for k in my calculation of t r . 

Rotational relaxation times, obtained from experiments are generally 0(10 sec). 

Therefore, it becomes convenient to express relaxation times on a different time scale. 

This time scale is determined by r, defined to be the mean time between collisions and 
-10 c 

is a number 0(10 sec). Relaxation times expressed on this time scale are known as 
collision numbers Z. Since we have two rotational relaxation times, t r and t r , we 
also have two rotational relaxation collision numbers 


r R 

R°T " - 


t 

ROT = ~ 




> 


(3) 


In a later portion of this report we will show that the calculation of Z r q t is essentially 
the calculation of the average change in rotational energy while the calculation of Z rQt 
reduces to the calculations of the average of the square of the rotational energy change. 
The following three sections of this report are devoted to this task. 
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Rotational Energy Change in a Collision 


The change in rotational energy produced by a collision will be calculated for the 
two-dimensional model of molecules constrained to the x-y plane. The molecules will 
have a mass M^., moment of inertia I k about an axis perpendicular to the x-y plane, 
and a permanent electric dipole moment // k constrained to the x-y plane. The inter- 
action will be the electric dipole -dipole interaction to which will be added a spherical 
hard core of diameter a. The situation is depicted in figure 2, which also defines the 
coordinates of the problem in the center of mass coordinate system. The Hamiltonian 
for the model is a sum of four terms, the relative translational kinetic energy H^, the 
rotational kinetic energy H^, a spherically symmetric short range interaction <f>, and 
the dipole -dipole interaction potential V. Using upper case letters for the translational 
coordinates and momenta, and lower case letters for rotational coordinates and mo- 
menta, we can write 


H = 


H rt + V 


(4) 


where 


h o =h t + h r + * 


(5) 


For simplicity, $ is taken to be a spherical hard core potential of diameter 
a = (a^ + a 2 )/2 where o ^ and are molecular diameters of the colliding mole- 
cules. The explicit expressions for Hrp, H^, and V are given in equations (6) to (8). 


H 


T 


1 

2m 



( 6 ) 


The coordinates and Q 2 are the radius and angle of polar coordinates, while 
and Pg are the conjugate momenta. 




(7) 


V = - [cos(q 1 - q 2 ) - 3 cosfaj - Q 2 )cos(q 2 - Q 2 )] 


( 8 ) 


(Qi)‘ 


8 



y 



Figure 2. - Unperturbed trajectories for scattering of point dipoies 
with impact parameters greater than molecular diameter. 


The coordinates and q 2 are angles specifying the orientation of the dipoles and p^ 
and p 2 are the conjugate momenta. All angles are positive when measured counter- 
clockwise and m is the reduced mass. The dipole -dipole interaction can be rewritten 
in the alternate form 


V = - — j^cos(q 1 - q 2 ) + 3 cos(q 1 + q 2 - 2Q 2 F) (9) 

2(Q 1 ) 3 

For the unperturbed problem, both H,p and are constants of the motion 

H° = E (10) 

H R = E R (11) 

where E is the relative translational energy, and the rotational energy. The ro- 
tational momenta are also constants of the motion for the unperturbed problem 
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I 


( 12 ) 


p£(t) = pj(f) 


while the rotational coordinates are given by 


q£(t) = qj(t’) + 



(t - 1* ) 


(13) 


The unperturbed problem for the translation part is also easily solved; however, for the 
purpose of this paper, we only need the time dependence of Qj and Q 2 . The time 
origin is tahen at the point of closest approach on the unperturbed trajectory (Q 2 = 0). 
For impact parameters b ^ a the incident particle is moving parallel to the y-axis and 
Q° and Q 2 satisfy 


r Q?(t)i 2 = — t 2 + b 

L J m 




cos Q 2 (t) = — - — b a a 

Qj(t) 




sinQ°(t ,= M\ 1/2 _1 
\ m / o°i 


Qj(t) 


(14) 


In equation (14), 2E/m is simply the square of the relative velocity. From these equa- 
tions it is apparent that Q° is an even function of t while Q° is an odd function of t. 
For impact parameters b < a, the unperturbed trajectories are again selected so that 
Q° and Q 2 have this same symmetry. This can be accomplished by a rotation of the 
coordinates x, y and the new situation is shown in figure 3. The Q° and Q° are de- 
scribed by the equations 
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Figure 3. - Unperturbed trajectories for scattering of point dipoles 
with impact parameters less than molecular diameter. 


In equation (15), the angle \p is defined by 

sin ip = sgn(t) (16) 

where sgn(t) represents the sign of the time variable. 

The change in rotational energy, AE^ = Hj^(t =+«>)- H^(t = - °°), that occurs as a 
result of the scattering can be calculated from equations (A18) and (A19) of appendix A 
by using H R as the dynamical variable F and by passing to the limit t — °° and 
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t 0 - 
since H 


R 


Since H R (t = -°o) = H R (t = -«>) by virtue of choice of initial conditions and 
is a constant of the motion along the unperturbed trajectory, we have 


AE r = + H^ 2) + . . 


QO 

-E 

n=l 


H ( n ) 

R 


(17) 


where the first and second order terms are given by the integrals 


H^P = lim 


f df n x (f)H°(t) 


t 0 - -oo 


H 


( 2 )_ 

R " 


lim 

t — CO 

to-- 


A/ 




> 


dt" 0 1 (t«»)n 1 (t>)H°(t) 


-J 


(18) 


The operator 0 X (t) is the Poisson bracket operator containing V. The dynamical vari- 
able H r contains only the rotational momenta P^ and not the translational momenta 
P^. Therefore, it is only necessary to evaluate the effect of the operator fi X (t') on the 
P, , and in turn, this requires the evaluation of the partial derivatives of equations (12) 

^ O 

and (13) with respect to p^ (t') 


Spjfo ) 

3p°(t») 


= 5 


kj 


aqg(t) 

3p°(t’) 




(19) 


and the partial derivatives of H R (t) with respect to p£(t) 


2 tt o i 


3H^(t) p£(t) 3 H R (t) 


_ 6 kj 


3pJ(t) 


J k 5p°(t)3p°(t) \ 


( 20 ) 


In equations (19) and (20), 5^ is the Kronecker delta. 

The calculation of rotational collision numbers Z r q t requires not AE R but AE R 
averaged over all possible initial conditions, or equivalently, averaged over all possible 
unperturbed trajectories. The unperturbed trajectories are determined by specifying 
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the rotational momenta p°, the rotational coordinates at t = 0, q£(0), the translational 

energy E, and the impact parameter b. For a given impact parameter, I shall use a 

Boltzmann type distribution function characterized by a translational temperature T and 

a rotational temperature T„. The normalized distribution function f can be written 

K 


f dE dp? dp? dq?(0) dq?(0) = ! 

1 i 1 i r(s-fl) 





e - H R/ k T R 

2 " kT R V¥2 



dq?(0) dq°(0) 

d Po — — 

2ir 2 7T 


0 < E < -co<p^<oo -7r < q^(0) < jt (21) 


The choice of the parameter s determines the type of distribution function, and it can be 
used to introduce an element of three-dimensionality into the calculation. Thus, s = 0 
corresponds to a two-dimensional Maxwell -Boltzmann distribution; s = 1/2 can be re- 
garded as either a three-dimensional Maxwell -Boltzmann energy distribution or a two- 
dimensional flux distribution; and the value s = 1 corresponds to a three-dimensional 
flux distribution. The distribution function (21) can be written in a slightly more con- 
venient form with the transformation 


2 \ V Z *1 


( 22 ) 


In terms of A^ the two -temperature distribution function takes the form 


f dE dA« dq?(0) dq°(0) = _L__ (IL)* e " E / kT &(*-) 
1 * r(s + 1) VkT/ VkT/ 


exp 




kT 


R 


(+> d<1 l (0) 


ffkT 


dA v ' dA 


R 


27 t 2ir 


(23) 


The distribution function, in the form (23), will be used to average (17). When 

AE d is averaged over the dipole orientation q°(0) and ql°^(0), all odd order terms in 
K i « th 

the expression (17) disappear. This is most easily shown by noting that the n order 

term will always contain a factor of the form 
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cos k q°(0) dqj(0) 



q°(0) dq°(0) 


k + i = n 


For the odd order terms, either k or i must be odd. But the result is always zero 
when an odd power of a trigonometric function is integrated over its period. Thus, if 
{ AE R ) denotes the average of AE R over the distribution function (23), we have 



( 24 ) 


The calculation of the rotational collision number Z r qt ret l uires 1116 calculation 
of ((AE r ) 2 ) . The quantity (AE r ) 2 can be expressed in terms of the perturbation ex- 
pansion (17) to give the equation 


«AEr) 


! )=EE <■«>> 

k=l z = 1 


(25) 


In this equation, as well as the equation for < AE R ) , all odd order terms vanish. Thus 
we get an expression for ((AE R ) 2 ), that is analogous to (24). 


«ae r ) 2 } = <(H^) 2 > + <(H^>) 2 + 2 hWh®). 


(26) 


If the first nonvanishing terms in (24) and (26) can be calculated, then we will have re- 
sults that are valid through third order in perturbation theory for both Z R Q,p and 
Z r qt- With these results we will be able to make a meaningful comparison of 
and Z R Qrp. That is, at least through third order in perturbation, we shall be able to 
make some conclusions regarding the equivalence or inequivalence of the two definitions 
of relaxation times t r , and r R . The next section will describe the calculation of 
while the following section will give the details in the calculation of <^(H^) 2 ^> 


Evaluation of Second Order Term 


The integrand of can be evaluated by the application of equations (20) and (A21). 

This gives the result 
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(27) 


o 1 (p»)n 1 (t')H°(t) = y [p^(t)s2 1 (t>')n 1 (t’)Pk(t) + 


It should be noted that the translational momenta P^ do not appear in (27) because 
is only a function of the rotational momenta P^. Furthermore, since the p£ and P^ 
are independent, (A22), when combined with equation (19), gives 


(28) 

<<t') 

Operating on equation (28) with f2^(t”) (or by using (A23) and remembering that the 
second partial derivatives of p£ vanish) gives the expression 


nVxiV)P?<t>.V 8v(t "> a2v<|,) (*' ~ *"> 

^ / rv n . o T 


/ ; aq°(t”) 9q°(t’)3q£(t’) 


X 3 


(29) 


Combining equations (28) and (29) with (27) and performing some rearrangements pro- 
duces 


^(f^O^H^t) 



im +(t - 

aq£(t’) 


3 2 V(t’) 

/ j h 9q°(t’)9q£(t’) 


(30) 


where the time label on p? has been suppressed since p? is a constant along the un- 
perturbed trajectories. 

Upon evaluating the partial derivatives of the dipole -dipole interaction potential, 
using the form (9), and averaging equation (30) over the orientations q°(0) and q°(0), 
there results the relatively simple expression 


*(<>>, E,b,t,t>) / dq°< 0 ) / dq^OJnV'JnVjH^t) = 

w 2 4 4 


fry* 
. 8I 1 : 2 . 


jcos[w ( '>(t" - t')] +9 cos p)(f. -t0 -2Q°(t'') + 2Q°(t’j] - t')sin[w ( ''(t'' -t’)] -9w (+) (t” -t')sin [w (+) (t" -t’) -2Q°(t") + 2Q°(t')] • 

3 

[Q°(t”)Q2(t’)] 


(31) 
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In this equation, I have used the notation 


a> (±) = ± P2^2 (32) 

I (±) = ^±12 (33) 

Equation (31) can be written in a somewhat more convenient form by introducing a para- 
meter e that will eventually be set equal to unity. Observing that 

(1 + 3/de) cos ecotl = cos cot - cot sin cot (34) 

|e=l 


we can define a function &(e) as 


#(e) = 


“2 2 t (+T 


8I i r 2 


(l + 3/3e) x 


cos^co( _) (t M -t’)]+ 9cos[^co( +) (t” -t’) -2Q°(t”) + 2Q2(t’)] 

|^Q j(t”)Q j(t’)j 3 


(35) 


such that #(1) = The integration of #(e) over the rotational energy distribution can 
be performed by using the relation between co^ and 

co (±) = ± (i/yfh ± * ( x Mh T i/v^ (_) ( 36 ) 

to obtain the result 
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To complete the evaluation of the second order term is still necessary to 

integrate equation (37) over the translational energy distribution and to perform integra- 
tions over t’ and t ,f between the limits given in (18). The integration limits for both 
time variables can be extended from to +°°, and the result multiplied by 1/2, be- 
cause of the symmetry of equation (37) in t* and t M . This then allows ^Ijj^to be 
written as 



lim 
e— 1 


2 2 ) 
MiM 2 

\4cr 6 kTy 



(1 + 3/3e) [5A(e) - 9B(e)] 


where the dimensionless functions A(e) and B(e) are defined by 


(38) 



(39) 



In equations (39) and (40), x = E/kT and the time scaling parameter Tj is defined by 
the expression 


T 


2 

1 


2I l l2 


I^kT 


R 


(41) 


Based on the equations for the unperturbed trajectories (14) and (15), two other time 
scales can be introduced by the definition 


2 

T 



2E x 


(42) 
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where 



,2 

rna 

2kT 


(43) 


The evaluation of the remaining integrals (39) and (40) is considerably simplified by 
the use of the Fourier transform. This is discussed in appendix B and the principal 
result is given as equation (B9). Combining equation (B9) with equation (38) gives 
the form 



2 2 

lim ^ 1/2 r(s+l/2) ^2 
e "" 1 r (s+l) CT 6 kT 



(44) 


Note that in equation (44) k is merely used as an integration variable and is not the 
volume viscosity. This is true for the balance of this section and also in the following 
section. The actual form of the density function p(k) can be obtained by evaluating the 
appropriate Fourier sine and cosine transforms. By performing the indicated operations 
on the parameter e, equation (44) becomes 


where 



ir 1/2 r(s+3/2) jfrf 
’ 2r (s + i) At 



pU)* 2 C 

r 2 2""I s+ ^/2 

Li + K 2 rJ 


c 


2 


T 1 _ / I 1 I 2 \ /T\ 
4t 2 " V^ma^ V T W 


(45) 


(46) 


The density function p(/c) depends parametrically only on b/a. The effect of the 
other parameters is given explicitly in equation (45). By using the letters C and S 
for the Fourier cosine and sine transforms, respectively, it becomes possible to write 
the density function p(/c) in the form 
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( 47 ) 


p(«) = 5C 2 [ — 18/ C 


sin 2 Q, 


f— ) cos 2 Qc 



sin Q 2 cos Q 2 


This expression can be evaluated analytically for b/a ^ 1, that is, when and Q 2 
are given by equation (14); however, for b/a < 1, p(/c) must be calculated numerically. 

In the former case (b/a > 1) all transforms are reducible to a form that can be recognized 
as either an integral representation of the modified Bessel functions of the second kind 
or a derivative of K^. The integral representation of is (ref. 25) 



dz 


cos k z 


i) oo ^ 1/2 

(z 2 + b/a 2 ) 


2 1/72 f—\ V K 

r(v+ 1/2) \2b/ "Iff/ 


(48) 


by using equation (48), p(/c) can be written as 


p(k) = 


/2\/a\ 4 

\ 7T/ \b/ 



2z 2 )z 2 K 2 (z) 


- 6z 3 K 1 (z)K 2 (z) + 2z 4 K 2 (z) 


- > 1 (49) 

a 


where z = /<b/a. This density function is finite everywhere and essentially goes to zero 
exponentially for large values of k. For impact parameters less than a the Fourier 
transforms that appear in equation (47) were evaluated numerically by the method de- 
vised by Longman (refs. 26 and 27). The function p(«r) is shown in figure 4 for several 
values of b/a. The value p(0) can be calculated from 



Figure 4. - Density function for rotational -translational 
energy transfer. 
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I 


2 

where A 
of b/a. 



( 50 ) 


1 - (b/cr) . The function /c p(fc) is plotted in figure 5 for several values 



Two integrations still remain undone. One of these is the average over impact pa- 
rameters and the other is the integration indicated in equation (45). Equation (45) can be 
evaluated analytically for impact parameters greater than a while for impact parameters 
less than a it must be calculated numerically. Designating the integral in equation (45) 

by £ , changing the integration variable to z = Kb /a, and substituting equation (49) for 
s 

p(k) gives 



20 



This integral can be evaluated by proceeding in a somewhat indirect fashion. The 
denominator is replaced by its Gaussian representation (ref. 28) by using 


/ oo 

d a „V2-1 


-cur 


This permits equation (51) to be written as 


•*s = -( £ 
TT lb 


— f 

r (s+3/2) / 

J 0 


/ 


v 2 2 


da a S+1/2 e~ a / dz e'"^^ z 
0 ~0 

Xz 2 |(5 + 2z 2 )z 2 K 2 (z) - 6z 3 K 1 (z)K 2 (z) + 2z 4 K 2 (z)] 


(52) 


The integration over z is of the type studied by Ragab (ref. 29) and can be expressed in 
terms of MacRobert’s E function (ref. 30) 


l 

77 - 1/2 


dz z* -1 K m (z)K n (z)e' z 


/i + m + n l + m - n £ - m + n I - m - n. I I + 1 . 


2 2 


: j3\ R J/3) > 0, R p (z ± m ± n) > 0 


By using this result, the remaining integral over a is essentially the Laplace transform 
of MacRobert’s E function. Transforms of this type were evaluated by Ragab (ref. 31) 
and have been tabulated by Erdelyi (ref. 32). Ragab showed that 

/ OO 

dX X k_1 e _X E^p; a r : q;p g : ^=E(p+l;a r : q;p g : z) 


where s k. Combining these two results with equation (52) gives for the integral jf g 
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where the E functions have been reduced to their lowest order. This expression can 
now be used in equation (45) to give 



For impact parameters less than a, £ and, therefore, the result corresponding to 
equation (54) must be obtained by numerical integration. The typical behavior of 
is illustrated in figure 6. It should be noted that, although ^l^^is continuous, it is 
sharply peaked at b/cr = 1. 



Figure 6. - Average change in dimensionless rotational 
energy as function of impact parameter. Mass distri- 
bution parameter, Z~ = 0.1. 
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Evaluation of the Second Order Term 



The evaluations of closely parallels the discussion of the previous section 

and, therefore, some of the details will be omitted. The integrand of can be ex- 

pressed in an appropriate form by combining the equations (A20), (A22), (19), and (20) 
to give 


oV)H°(t) = - 



P°(t) 


8V(t t ) 


*k 3q°(t') 


k=l 


ah 


This may now be used to express (H^O in the form 


(55) 



2 

) 





dt „ Pk p k av(t") 9V(t > ) 

r k 8q°(t”) 3q°(t’) 


(56) 


Upon evaluating the partial derivatives of the dipole -dipole interaction potential (9) and 
averaging the integrand of (H^)^ over the orientations q°(0) and qg(0) we obtain the 
function 
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■ ■nil 


i 


This function although similar to the function defined in (31), is not identical 
to The next step in the calculation is to average H over the rotational energy dis- 
tribution and we then obtain 



x 


•fl + 9 cos 



2 (t’) -Q 2 (t* 



exp 


e 2 (t” - t') 2 I^kT R 


[Ql(t”)Qi(t’)] 3 


(58) 


The right side of this expression is identical to the right side of equation (37) apart from 

a factor 2kT„. If we now combine equation (58) with the results of the previous section 
K 

we obtain the expression 


^Hr 1 ’) 2 ) - 2kT R<?R 2) >“ 


_ 7T 1/2 r( s+ 3/2) M 


r (s+ 1) 


/ oo 

dK pW! 

(' ♦ 


£ 

>\ s+3/2 

1 


(59) 


The Calculation of Z R qj and Z R qj 


Equation (1) can be used as a basis for the calculation of r R . The derivative 
dEra/dt represents the change in rotational energy per collision averaged over the num- 

XV 

ber of collisions per unit time. This is equivalent to the average change in rotational 

energy per collision, /^E R \ divided by r , the mean time between collisions. Thus, 

\ *vAv c 
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I 


at t = 0 equation (1) becomes 



E r (T) - E R (t = 0) 
t R 


( 60 ) 


Assuming the molecules to be rotationally unexcited at t = 0 and solving for t c t r ^ we 
find 


-1 

ROT " 


t c t R 


1 

Er(T) 



( 61 ) 


The calculation of /^EjX includes only the contribution of trajectories with impact 

parameters out to some maximum impact parameter b. All of the impact parameter 
dependence is contained in the integral £ defined in equation (51). Because £ goes 
to zero so rapidly, the integration can actually be extended to infinity with little error. 
We, therefore, define a function 0 by 

toil 



where n can either take the value 0 or 1. The value n = 0 corresponds to a two- 

dimensional calculation while n = 1 corresponds to the M three -dimensional" case. In 

terms of 0 , we can now write the expression for the change in rotational energy 

sn 

averaged over orientations, translational and rotational distribution functions, and impact 
parameters in the form 
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TABLE II. - NUMERICAL VALUES FOR 


FUNCTION O a „ 
sn 


Mass 

distribution 

parameter, 

? 2 

Average rotational energy 
change per collision 

e (l/ 2 )o/ 0 ll 

0 (l/2)O 

©ii 

0 

0 

0 

1. 4443 

2. OxlO” 3 

2. 9884X10” 1 

2. 0691X10" 1 

1.4207 

4. OxlO' 3 

3.8473 

2.7081 

1.4207 

6. OxlO' 3 

4.3908 

3. 1215 

1.4066 

8. OxlO” 3 

4.7823 

3.4228 

1.3972 

1. OxlO” 2 

5. 0827 

3.6554 

1. 3905 

2. OxlO' 2 

5.9327 

4.3126 

1. 3757 

3. OxlO” 2 

6.3115 

4. 5954 

1.3734 

4. OxlO' 2 

6.4967 

4. 7225 

1.3757 

5. OxlO” 2 

6. 5819 

4. 7698 

1.3799 

7. 5X10" 2 

6. 5876 

4. 7268 

1.3937 

1. OxlO” 1 

6.4675 

4. 5916 

1.4086 

1. 5X10' 1 

6. 1259 

4.2626 

1.4371 

2. OxlO" 1 

5.7738 

3.9464 

1.4631 

2. 5X10' 1 

5.4490 

3. 6659 

1.4864 


Table II gives the function Q an in tabular form. 

oil 

To assist in the discussion of the function 0„ , let us reexamine equation (45). The 
P sn 

parameter q that appears in this equation is proportional to the parameter used by 

Sather and Dahler (refs. 13 and 14) in their rough sphere calculations. It differs from 

o 

their parameter 41/mcr only by an additional factor of (T/4T r ). Therefore, the inte- 
grand of equation (45) has a moment of inertia dependence quite similar to the rough 
sphere results. This is particularly true for the case s = 1/2. Because of this fact and 
because of the sharply peaked behavior of £ , it is desirable to define a comparison func- 

tion 0° , analogous to 0 , but based on the use of a density function k p (k) = 

sn. 7 sn. 

5{k - 1)6(1 - b/o) where 6(/c) is the Dirac delta function. This gives 


o _ (n + IK 
sn 


9 s+3/2 

(i + n 


The function 0° can be thought of as corresponding to the rigid sphere result for the 
sn 


sn 


model used in this paper. For comparison purposes, figure 7 gives plots of both 0 

and 0° . From figure 7, it can be seen that averaging over the density function p(x) 
sn 2 

has the effect of shifting the maximum from relatively large values of £ in the function 
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.7 



0° to relatively small values of If 2 in the function ©„ . The net result is that 

sn 

Z ROT’ for this model » can be either an increasing or decreasing function of I depend- 
ing upon whether £ 2 lies to the left or to the right of the maximum. 

We are now in a position to calculate z r qt from equation (61). The mean time 
between collisions, r , is given by 



(3 -dimensions) 


(2 -dimensions) 


(64) 


where v is the number density with dimensions appropriate to either two or there dimen- 
sions. Substituting this equation, together with equation (63) into equation (61) gives 
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/dE 


R\ 


Z ROT " t c t R “ T c 


dt 


t=0 


E r ( T ) 


w 1//2 r(s + 3/2 ) 
2r(s+ 1) 


2 2 


6 


a kTE R (T) 



(65) 


In the usual case the two colliding molecules are identical and, therefore, \i ^ = jug = M 
and a 2 = a * Also, we can write E^ = ykT/2 where y is the number of rotational 
degrees of freedom and is either 2 or 3. Thus, equation (65) can be written as 


-1 _ 1 1/21(8+3/2) 

ROT V r(s+ 1) 



9 S „(? 2 ) 


( 66 ) 


The adjustable parameter b can be eliminated by considering a ratio of collision num- 
bers. If, in addition, we consider the near equilibrium situation, we may set T/T^ = 1 
to obtain 

Z ROT^ 1, T ) _ y _i / M 2\ £l\ /T \ 2 0 sn^2^ ( 67 

Z ROT^ 2 ’ T o^ r 2 V*l/ \ 2/ W © sn (^) 


where the notation z r ot^’ T ) re P resents the collision number of species 1 at a temper- 
ature T. 

The equation (2) will serve as a basis for calculating z r qt* Tlie classical equiva- 
lent of the Wang Chang, Uhlenbeck, and DeBoer expression (ref. 24) for the volume 
viscosity k can be written in our notation as 



Combining this equation with equation (2), the three-dimensional expression for t q (65), 
and using yk/2 for C.^ gives 
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2 


(69) 


(rv' N-l _ (n + 1) 
< Z ROT> - — ' " 



n = 1, s = 1, T/T r = 1 

Strictly speaking, the parameters n and s were unnecessarily introduced into this 
equation since the Wang Chang, Uhlenbeek, and DeBoer formula for k requires both of 
these parameters to be unity. However, if we relax this restriction and permit the com- 
bination n = 0, s = 1/2 we will have a greater parallelism between the calculation of 
? 1 
Z ROT and Z ROT* The collision number Z r q t can then be given an analogous inter- 
pretation. That is, since ^[AE R ) 2 ^was estimated for a two-dimensional system, n = 1, 
s = 1, corresponds to a quasi three-dimensional calculation while n = 0, s = 1/2 repre- 
sents a two-dimensional calculation. 

We can now substitute for ^AE R /kT R ) 2 ^> in equation (69) by using its perturbation 
approximation given in equation (59). If we then combine this result with the definition of 
the function 0 ori given in equation (62) we find 

oil 


( Z ROt) 


1 7T 1//2 r(s+3/2) ^1^2 q 
2y r(s+l) CT 6 (kT) 2 sn 


(70) 


Setting n 1 = p 2 and comparing this with the expression for Z RQT given in equation (66) 
shows that within the level of approximation employed in the calculation we obtain the 
surprising result 


Z ROT ~ 



(71) 


Thus, both def ini tions give identical predictions about the dependence of the rotational 
collision number on the temperature, the dipole moment, and the moment of inertia. 


COMPARISON WITH EXPERIMENT 

It would be unrealistic to expect equations (66) or (70) to agree quantitatively with 
experiment in view of the approximate nature of the calculation. The best that one can 
hope for is a relative comparison based on equation (67), which holds equally well for 



Z ROT and Z ROT t> ecause °f equation (71). Hence, attention will be confined to a com- 
parison of equation (67) with experimental data. 

The molecular diameter and the moment of inertia must be evaluated before the 
theoretical results (eq. (67)) can be compared with experimental data. The moment of 
inertia will be taken to be the arithmetic mean of a molecule’s principle moments of 
inertia. O’Neal and Brokaw (ref. 33) suggest that a be set equal to the viscosity colli- 
sion diameter cr ■ and calculated by the relation 



5 V?rMkT 
16 7T77(T) 


(72) 


In expression (72), rj( T) is the viscosity at a temperature T. O'Neal and Brokaw used 
a only in the calculation of t ; I shall use the identification (72) for a wherever it 
appears in the formulae. This use of a for <r might be viewed as an attempt to account 
for the fact that the spherically symmetric short range repulsion, although rising 
abruptly, does not rise vertically as in a rigid sphere. 

By using equation (72) and considering first the quasi-three-dimensional result 
(s = 1, n = 1), equation (67) takes the form 


z ROT< 1 » T > n M 4 /S\ /T\ 3 p2< T o>] 2 e il (2 ’ T 0> 

Z ROT< 2 ’ T o> y 2 \^2/ \ M 2/'V k (T) J e il (1 ' T > 

2 

The variable % , required for the evaluation of 0jp is to be calculated as 




2 


16tt 



5 VjrMkT 


(73) 


(74) 


Some qualitative comments can now be made concerning the behavior of equation (73). 


According to this equation, an increase in dipole moment decreases the collision number. 

o 

Furthermore, since all the moment of inertia dependence is contained in £ and since 
2 

the values of £ (as calculated from eq. (74)) lie to the left of the maximum in 0^, an 


increase in I causes a decrease in 


Both of these results are in accord with 


experimental observations. Finally, Z increases with temperature when calculated 
by equation (73), however, because of the compensation due to viscosity, the increase is 
at a considerably lower rate than indicated by the explicit T 3 dependence. Based on 


some typical calculations, the increase is very nearly linear at the lower temperatures 
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but becomes somewhat steeper as the temperature is increased. At present, there is 
sufficient uncertainty about the temperature behavior of z rqt so that one cannot con- 
clude that equation (73) either agrees or disagrees with experiment. 

A comparison between equation (73) and experimental results is made in tables III 
and IV. To evaluate equation (73), it was necessary to use the data of table I and, in 
addition, viscosity data. The viscosity data cited by Baker and Brokaw (refs. 1 and 2) 
were used in the evaluation of equation (73) for all molecules except SOg, HgS, HF, DF, 
and the fluoromethanes. For SOg and HgS the viscosity was obtained by fairing a curve 
through the data contained in the references cited by Svehla (ref. 34). The HF viscosity 
was calculated by using the extrapolation formula of Posey (ref. 35). The viscosity cross 
section of DF was taken to be equal to HF. The viscosity data for the fluoromethanes was 
taken from Svehla’ s tabulation (ref. 34). 

The theoretical temperature dependence of the rotational collision number for SO 2 
and NHg are compared with experimental values in table HI. Only these molecules could 


TABLE III. - COMPARISON OF THEORETICAL 
AND EXPERIMENTAL TEMPERATURE 


DEPENDENCE OF Z RQT 


Tempera- 

ture, 

K 

Z(T)/Z(T q ) 

Theoretical 

Experimental 

Two- 

dimensional 

Three- 

dimensional 


so 2 


300 

1. 00 

1. 00 

p 

♦— k 

0 

0 

500 

1.42 

1. 61 

1. 1 to 2. 1 

700 

2. 00 

2. 42 

1. 9 to 3. 9 

900 

2. 73 

3. 42 

3. 7 to 9. 4 

NH 3 

311. 2 

1. 00 

1. 00 

0 

0 

1—1 

rQ 

333. 2 

i-* 

0 

CO 

1. 05 

1. 02 

353. 2 

1. 05 

1. 09 

1.03 

373. 2 

1. 07 

1. 13 

1. 12 

393.2 

1. 09 

1. 17 

1. 20 

413.2 

1. 12 

1. 21 

1. 25 

433. 2 

1. 14 

1. 25 

1.29 

453. 2 

1. 16 

1. 29 

1. 38 

473. 2 

1. 18 

1. 33 

1. 59 


a Ref. 4. 
b Ref. 5. 
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TABLE IV. - COMPARISON OF THEORETICAL AND EXPERIMENTAL 


Zrqt for A GIVEN TEMPERATURE 



Tempera- 


Theoretical 


Experimental 


ture, 



1 




K 

Two- 

Three- 



dimensional 

dimensional 


Z DC1 /Z HC1 

300. 1 

b 0. 

54 

c o. 

76 

a 0. 43 


328. 5 




77 

. 50 


375.8 





, 56 


423. 1 





, 60 


471.4 





, 62 

z d 2 o/ z h 2 ° 

381.2 
426. 1 

b 0. 

63 

c 0. 

75 

76 

d 0. 52 


478. 0 




75 



525.6 




76 

— 

Z NDo/ Z NH, 

300.0 

b 0.64 

c o. 

78 

a 0. 69 

o/ o 

329. 1 




79 

. 70 


374.6 





. 72 


424.2 





. 72 


424.2 





. 74 


474. 5 




.. 

. 77 

Zl5 NH 3 / ZNH 3 

300.0 

424.2 

b l. 05 
1. 05 

c l. 02 
1. 01 

__ 

a l. 05 
1. 04 

Z Df/ Z HF 

373.8 

b 0. 58 

c 0. 78 

e 0. 44 

422.3 

. 58 

. 79 

.41 

Z(H 2 S, 298. 3° K)/ 



c 14. 8 


f 

*11.6 

Z(H 2 0, 381. 2° K) 







Z(H 2 S, 298.3° K)/ 
Z(HC1, 300. 1° K) 



c 2 

.9 


g 5. 0 

Z(CH 2 F 2 , 300° K)/ 
Z(CHF 3 , 300° K) 



c 

. 56 


h . 77 

Z(CH 2 F 2 ,300° K)/ 



c l. 08 


CO 

1 

Z(CH 3 F, 300° K) 








a Ref. 2. 

^Modified Sather and Dahler (refs. 2, 13, and 14); 

Z ROT a («/MoJ)/(l + 4I/Ma^) 2 . 

C Eq. (60). 
d Ref. 1. 
e Ref. 3. 
f Refs. 1 and 6. 
g Refs. 2 and 6. 
h Ref. 10. 
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be compared with theoretical results because the experimental data for the other mole- 
cules indicate that decreases with temperature. To put table ni in a somewhat 

better perspective, it should be pointed out that the collision numbers obtained from the 
ammonia data of Baker and Brokaw (ref. 2) have a temperature dependence that is oppo- 
site to those of Srivastava and Das Gupta (ref. 5). Thus, it is somewhat surprising that 
the experimental and theoretical values agree as well as they do. In the light of the 
present uncertainties concerning the experimental data, this agreement must be viewed 
as largely fortuitous. 

A comparison of isotope and dipole moment effects, as calculated from equation (73), 
is given in table IV. In addition, the results from the modified Sather and Dahler formu- 
la (ref. 33) are also given. Here, the agreement is not as good as was the case for 
table IV, in fact, some of the calculated numbers differ from the experimental results 
almost by a factor of two. In general, it can be said that equation (73) qualitatively 
mirrors the isotope and dipole moment effects, however, quantitative agreement between 
theory and experiment is lacking. 

In addition to choosing values of a and I, it is necessary to decide what values 
the indices s and n should take. The qua si -three -dimensional case corresponds to 
the choice s = 1 and n = 1 while the two-dimensional case corresponds to the choice 
s = 1/2 and n = 0. Up to now, we have only considered the ’’three-dimensional case." 
The relation of the two choices can be easily established from equation (67). 


Z ROT^ 1, T ) 

_ Z ROT^ 2 ’ V 
Z ROT^ 1 ’ T ) 

Z ROT^ 2, T o) 

2 -2 

The parameter £ is generally less than 5x10 and, thus from table II, it follows that 
the 0 gn in equation (75) produce differences of the order of 1 percent. The ratio of 
viscosity cross sections, on the other hand, can produce differences of the order of 
10 percent or more. Since the viscosity cross sections of the isotopic pairs are gener- 
ally taken to be equal at a given temperature, we can conclude that both the two- and 
’’three-’* dimensional cases will predict essentially the same moment of inertia depend- 
ence. However, because of the rather strong temperature dependence of the viscosity 
cross section, there will be a substantial difference in the predicted temperature de- 
pendence. In general, the two dimensional case displays a considerably slower increase 
of z rot witil temperature. This can be seen from the data in table HI. It does appear 
that the "three” dimensional case agrees with the experimental results somewhat better 
than does the two-dimensional case. 


[ 2 -dim _ _1 

On 


© 


,( 2 ) 


1/2 O' 


© n (2) 


e il^) 

0 l/2 0<» 


(75) 


3 -dim 
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CONCLUDING REMARKS 


The calculations that were described in the previous sections show that the two 
definitions of a relaxation time are equivalent for our model, at least through third order 
in perturbation. Furthermore, a comparison with experimental values given in the liter- 
ature seems to indicate that the calculated rotational collision number is in moderate 


agreement with experimental results. 

The reason for the lack of better agreement between the theoretically and experi- 
mentally deduced effects of temperature, dipole moment, and moments of inertia on 
Zrqt is problematical. In part, it may be due to the fact that the theoretical calcula- 
tion was essentially two dimensional and not three dimensional as required by reality. In 
recent Monte Carlo calculations both Bunker and Blais (ref. 36) and Karplus and Raff 
(ref. 37) detected considerable differences between two- and three-dimensional results. 

A portion of the disagreement might possibly be a result of the data used in evaluating 
equation (73). For example, a 2 percent error in the viscosity ratio results in a 4 per- 


cent error in the ratio of 
an 8 percent error in the 


Z ROT* 

Z ROT 


A similar error in the dipole moment ratio produces 
ratio. Such errors are undoubtedly contributing factors. 


For example, the viscosity cross sections of isotopic pairs were generally taken to be 


equal. Only in the case of the pair HgO and DgO was any correction made. A third rea- 
son for the discrepancies is the apparent sensitivity of Z r q t to relatively small exper- 
imental errors in thermal conductivity. Finally, some of the difficulties might possibly 


be a result of some deficiency in the theory used for extracting Z r q t from thermal 
conductivity measurements. At present, it is impossible to assess the relative impor- 


tance of each of these factors. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, April 18, 1969, 
129-01-06-01-22. 
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APPENDIX A 


CLASSICAL PERTURBATION THEORY 

Garrido (refs. 21 and 22) has presented a formulation of classical mechanics in the 
interaction picture. Formally, the development is very similar to the corresponding 
treatment of the Dirac or interaction picture of quantum mechanics. For any classical 
system with the Hamiltonian 


H=H 0 + V 


(Al) 


the solution to Hamilton's equations can be written in the form 


q = q(p o> q o> - V 

P - 0(Pq, qQ, t - tg) 


(A2) 


where q^ and Pq are the values of the coordinates and momenta at the initial time t^. 
The solution of Hamilton's equations for the unperturbed systems, characterized by the 
Hamiltonian Hq, can be written in a similar form 


q = q°(p 0 , q 0 > 
p = p°(p 0 , q Q , 


> 


t-t 0 ) 


(A3) 


Thus, it is assumed that the initial values of the coordinates and momenta for both the 
perturbed and the unperturbed system are identical. 

The time evolution of any dynamical variable 


F = F(q, p, t) 


can easily be obtained for either the perturbed or the unperturbed system by substitu- 
tions of the appropriate solutions (A2) or (A3). For the perturbed system, 

F(t) = F jq(p 0 , q Q , t - t Q ), p(p Q , q Q , t - t Q ), t] s U(t, t Q )F(q 0 , p Q , t Q ) (A4) 

and, similarly, for the unperturbed system, 

F°(t) = F [q°(p Q , q Q , t - t Q ), p°( P() , q Q , t - t Q ), t] = U°(t, t Q )F(q 0 , p Q , t Q ) (A5) 
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The evolution operators U and U° defined in equations (A4) and (A5) have the usual 
properties 


U(t, t")U(t", t‘) = U(t, t r ) ' 
U(t, t) = 1 > 

U _1 (t, V) = U(t*, t) 


(A6) 


From a consideration of infinitesimal time translations, one obtains the differential 
equations satisfied by the time evolution operator 

dU(t, t n ) 

—V- = 0(t)U(t, t Q ) = U(t, t 0 )J2(t 0 ) (A7) 

where the operator £2(t) is defined by 

n(t) = i--[H,] (as) 

9t 3t 9p 0q 0q 0p 


The second equality in equation (A7) follows from the fact that for any operator O(t) 

O(t) = U(t, t 0 )O(t Q )U(t 0 , t) (A9) 

Analogous equations can be written for the unperturbed system. The solution of equa- 
tion (A7) can be easily obtained from the integral form of this equation. 


u(t, t 0 ) = 1 + dPO(t')U(t», t 0 ) = 1 + f dt' U(V, t 0 )Q(t 0 ) (A10) 

t 0 


Using the first equality in equation (A10), we get as the iterative solution 


U(t, t Q ) = P exp 



(All) 


where P is the Dyson time ordering operator. Using the second equality in equation 
tion (A10), we get the alternate result 
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u(t, t Q ) = exp[(t - t 0 )^(t Q )] 

Therefore, the operator U(t, tQ) is nothing more than a Taylor expansion about tp. 

To apply these results to perturbation, it is only necessary to observe that equa- 
tions (A4) and (A5) can be combined to give 

F(t) = U(t, t 0 )F(q 0 , p 0 , t Q ) = U(t, t 0 )U°(t 0 , t)U°(t, t 0 )F(q 0 , p Q , t Q ) 


or 


F(t) = S(t)F°(t) (A12) 

where the operator S(t) is defined as 

S(t) = U(t, t 0 )U°(t 0 , t) (A13) 

From the differential equations for U and U° (eq. (A7)), one readily obtains the differ- 
ential equation satisfied by S 


^=S(t)f2 1 (t) (A14) 

dt 

where the operator f2*(t) is a Poisson bracket 

f2(t) - £2°(t) = oV) = - fv, ] = ^ — ' (A15) 

3p°(t) 3q°(t) 0q°(t) 3p°(t) 


and V is evaluated along the unperturbed trajectories. The iterative solution of the in- 
tegral equation obtained from equation (A14) can be written as 


S(t) = 9 exp 



dt» ftV) 


1 + 


+ 



dt’ fiV) 


dt’ 



dt” fi 1 (t")f2 1 (t T ) + . 


(A16) 
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The time ordering operator & in equation (A16) orders time -dependent operators in 
sense opposite to that of the Dyson time ordering operator P. 




= 


t’ < t” 


t” < t’ 


(A17) 


The relation between F(t) and F°(t) given by equation (A12) can now be rewritten by 
using the iterative expansion of the operator S(t) given by equation (A16). 


F(t) = 


z 


F^(t) = + F^t) + F^(t) + . . . 


(A18) 


where the first three terms are given explicitly by 

f(°V) = F°(t) 


/ 


F^(t) = / dt’ D 1 (t’)F°(t) 


L 0 


F< 2 >(t) 


l 0 l 0 


dt” n 1 (t”)n 1 (t’)F°(t) 


(A19) 




To evaluate the integrals given in equation (A19) one must be able to determine the 
effect of f2^(t’) on the function F°(t). This is possible when one knows the solution to 
the unperturbed problem in the form of equation (A3). Using equations (A5), (A15), and 
the chain rule for differentiation, we obtain 


= — - 


3p°(t) 


9F°(t) 3q°(t) + 0F°(t) 0p°(t) 


Sq°(t) 0q°(t>) 9p°(t) 0q°(t») 


av 

aq°(f) 


0F°(t) 0q°(t) ; 0F°(t) 0p°(t) 
3q°(t) 0p°(f) 0p°(t) 3p°(f) 
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or, upon rearranging, 


S2^(t’)F°(t) = dF (t) 3F (t) 


9q°(t) 


n (t’)q (t) + 


3p°(t) 


^(t’)p U (t) 


Explicitly indicating that more than one set of coordinates and momenta may enter the 
problem, we have for the special case where F depends upon momenta alone, the result 


S2 V)F°(t) =? 



SO&nWo, 


3p£(t) 


n (t’)pj(t) 


(A20) 


The second order term is somewhat more complicated 


n 1 (t ,, )n 1 (t , )F u (t) 


»)F°(t) X n 1 (t M )n 1 (t , )p?(t) 

/ ^pjfo) 



S2 1 (f * )p°(t)S2 1 (t* )p£(t) (A21) 


k,l 


3p°(t)Sp°(t) 


When the perturbation only depends on coordinates, we can write 


nV)p£(t) = - 



and the second order term 


o V)« V)$t) =\ 


3q°(t") 3p°(t») 



8V(t') Sp ifo 

0q°(t») 3p°(t») 


9V ( t t) dpjfo) 

3q°(f) 3p°(t») 


av(t») s 2 v(t t ) 8q ° (tt) ap k^ | \ 9V(t") av(t') 

1— i 3q?(t”) 9q°(f)9q°(f) 9p°(t") 9p>) / J 3q?(t») 3q°(t») 9p?(t»)3P>) 
4 7 i J 1 L J L A 1 A t' 


3 2 Pk(t) 




hi 


(A22) 


(A23) 
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APPENDIX B 


EVALUATION OF A(«) AND B(«) 

The function A(e), defined in equation (34), contains an integral of the convolution 
or Faltung type 





dt x dt 2 f(t 1 )g(t 2 )r(t 1 - t 2 ) 


(Bl) 


where r(t) is the Gaussian function 


r(t) = exp 



(B2) 


The function B(e), defined in equation (35), will also contain a linear combination of 
integrals of the type in equation (Bl) if the trigonometric function in its integrand is 
expanded. The convolution theorem (ref. 38) can be used to rewrite equation (Bl) in 
terms of the Fourier transforms of the functions in its integrand 


(2jt) 1/2 



d/cF*(K)G(/c)R(/c) 


(B3) 


where the asterisk is used to denote complex conjugation and where the Fourier trans- 
forms are defined as 


G(k) = (2 rr)- 1/2 



dt g(t)e iKt 


(B4) 


For the case at hand, the functions f(t), g(t), and r(t) possess a definite symmetry being 

either even or odd functions of t. Furthermore, f(t) and g(t) will always have the same 

symmetry. For an even function the Fourier transform is identical to the Fourier cosine 

transform while for odd functions it differs from the Fourier sine transform by a factor 
1/2 

of i = (-1) ' . The sine and cosine transforms are defined by the expression 
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(B 5) 



Thus, for example, the Fourier transform of the Gaussian (B2) is 


R(k) = R c (fc) = 



cos ict = 2 


1/2 


' Il \- ( * Tl/2e)2 

2e/ 


(B6) 


Combining this result with equation (B3) and using the symmetry properties of the sine 
and cosine transforms gives 



In the evaluation of A(e) and B(e) it is still necessary to integrate integrals of the 
type (eq. (B7)) over x = E/kT. This can be accomplished without knowing the functional 
form of the transforms F and G. This is possible because the functions f and g de- 
pend only on the combination (t/r) where r contains all the x dependence (eq. (37)). 
This means that with an appropriate change of integration variable equation (B5) takes 
the form 



Using this equation, together with an analogous expression for the transform F, and 
making the change of variable kt — k gives 



41 



Now, averaging y(x) over the energy will give 
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